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The  transfer  function  associated  with  the  Knox-Thompson  speckle  imaging  technique  is  investigated.  Numerical 
model  transfer  functions  using  log-normal  statistics  for  perturbations  of  the  complex  wave  front,  the  near-field 
approximation,  and  a  Kolmogorov  spectrum  for  atmospheric  turbulence  statistics  are  presented.  Simple  approxi¬ 
mations  for  the  transfer  function  are  discussed.  As  with  the  transfer  function  of  Labeyrie’s  speckle  interferometry 
technique,  the  portion  beyond  the  seeing  limit  can  be  represented  as  the  transfer  function  of  an  unaberrated 
telescope  times  a  seeing-dependent  constant.  An  additional  factor  depends  on  the  frequency  shift  of  the  Knox- 
Thompson  cross  spectra.  The  influence  of  the  frequency  shift  on  the  reconstructed  phase  error  is  discussed  for 
simple  reconstruction  problems. 


1.  INTRODUCTION 

The  Knox-Thompson1  (KT)  speckle  imaging  technique 
plays  a  significant  role  in  present  methods  of  high-resolution 
optical  astronomical  imaging.  Speckle  imaging  techniques 
permit  the  recovery  of  small-scale  information  on  astronom¬ 
ical  sources  that  normally  is  destroyed  by  the  deleterious 
effect  of  turbulence  in  the  Earth’s  atmosphere  on  the  imag¬ 
ing  process.  Unlike  Labeyrie’s2  speckle  interferometry, 
which  measures  the  autocorrelation  of  the  observed  intensi¬ 
fy  distribution,  the  KT  technique  permits  the  recovery  of 
the  full  information  on  the  observed  intensity  structure. 
Because  the  technique  is  relatively  simple  to  apply  com¬ 
pared  with  other  techniques,  such  as  speckle  masking,3  it  has 
found  wide  application  in  high-resolution  imaging  astrono¬ 
my,  ranging  from  stellar  sources4  to  extended  objects,  such 
as  the  Sun.3 

In  speckle  interferometry  and  speckle  imaging,  a  large 
number  of  short-exposure  pictures  of  the  source  of  interest 
are  taken.  Each  picture  is  affected  by  atmospheric  turbu¬ 
lence,  and  if  the  exposure  is  short  compared  with  the  atmo¬ 
spheric  correlation  time  (of  the  order  of  10  msec),  the  ob¬ 
served  intensity  distribution  of  the  ith  picture,  may  be 
represented  by  the  true,  unaberrated  intensity  distribution, 
/o,  convolved  by  a  randomly  speckled  point-spread  function, 
PSF„  which  typically  spreads  over  1  arcsec: 

/,(x)  =  70(x)  *  PSF,(x),  (1) 

where  x  denotes  a  two-dimensional  focal-plane  coordinate 
and  *  denotes  the  convolution  operation.  When  long  expo¬ 
sures  are  taken,  the  speckles  average  out,  and  the  effective 
point-spread  function  is  smooth  and  does  not  permit  the 
recovery  of  detail  smaller  than  its  average  diameter  xs,  called 
the  seeing  limit. 

The  speckles  in  the  point-spread  function  are  of  the  size  of 
the  theoretical  resolution  element  of  the  instrument,  and  a 
suitable  statistical  analysis  of  their  modulation  permits  the 
recovery  of  object  structure  information  up  to  the  resolution 
limit.  Such  an  analysis  is  frequently  done  in  the  Fourier 


domain.  If  F;  and  F0  are  the  two-dimensional  Fourier  tran 
forms  of  7,  and  7o,  respectively,  then 

F,(f)  =  F0(f)S,(f),  0 

for  which  S,  denotes  the  instantaneous  optical  transfer  fun 
tion,  i.e.,  the  Fourier  transform  of  PSF,(x),  which  depenc 
on  the  telescope  and  the  atmosphere,  and  f  denotes  a  two 
dimensional  spatial  frequency.  The  essence  of  the  KT  tech¬ 
nique  is  to  evaluate  the  average 


±2Fi(r>Ff(  f-A) 


1 

=  F0(f)F0*(f  -  A)  X  ~  £  S,(f)S,*(f  -  A)  (3) 


for  a  fixed  spatial-frequency  spacing  A.  N  is  the  number  of 
Fourier  transforms  processed.  Both  sides  of  Eq.  (3)  are 
complex  functions.  There  is  some  controversy  about  what 
to  call  quantities  such  as  F(f)F*(f  —  A);  the  term  cross 
spectrum  seems  to  prevail,  so  I  shall  use  it  throughout  this 
paper. 

Equation  (3)  shows  that  the  observed  average  cross  spec¬ 
trum  is  given  by  the  object  cross  spectrum  multiplied  by  an 
average  telescopic-atmospheric  term,  the  cross-spectrum 
transfer  function  (CTF).  Knox  and  Thompson1  realized 
that,  in  the  ensemble-average  limit  N  -*■  the  CTF  is  a 
purely  real  function  as  long  as  A  is  small  compared  with  the 
seeing  cutoff  frequency  fs,  which  is  reciprocal  to  the  seeing 
limit  x,.  The  Fourier  phase  of  the  observed  cross  spectrum 
is  equal  to  the  phase  of  the  object  cross  spectrum;  the  object 
Fourier  phase  may  be  recovered  from  the  observed  cross- 
spectrum  phase  by  means  of  suitable  integration  techniques. 
If  the  CTF  is  known,  Fourier  amplitudes  may  be  also  recov¬ 
ered  from  the  cross  spectra  [Eq.  (3)];  amplitudes  and  phases 
can  be  recombined  and  inversely  transformed  to  obtain  the 
reconstructed  object  intensity  distribution.  Because  of  the 
better  signal-to-noise  ratio,  Fourier  amplitudes  are  com¬ 
monly  obtained  by  using  Labeyrie’s2  method  of  analyzing 
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cross  spectra  for  A  =  0  (referred  to  as  the  Labeyrie  case  in 
this  paper). 

The  accuracy  to  which  phases  can  be  recovered  from  ob¬ 
served  cross  spectra  depends  on  A,  because  the  amplitudes 
of  the  observed  cross  spectra  decrease  as  A  increases.  As  a 
rule  of  thumb,  Knox  and  Thompson1  require  that  A  be 
smaller  than  0.5  /„•;  Fried8  requires  that  A  <  0.2  fs.  A  more 
accurate  description  of  how  cross  spectra  amplitudes  are 
affected  by  the  choice  of  A  is  required  in  order  to  study  the 
dependence  of  phase  errors  on  A.  Also,  the  choice  of  A 
imposes  a  lower  limit  on  the  size  of  the  observed  field  of  view. 
If  an  extended  source,  such  as  the  solar  surface,  is  observed, 
a  small  field  of  view  is  required  in  order  to  match  the  size  of 
the  isoplanatic  patch.6-7  Once  the  field  size  is  selected,  A 
cannot  be  made  smaller  than  the  fundamental  frequency 
given  by  that  size,  and  it  is  conceivable  that  under  unfavor¬ 
able  seeing  conditions  the  KT  algorithm  will  fail.  Also, 
variations  of  the  KT  technique  use  sets  of  cross  spectra  with 
varying  A,8  and  in  these  cases  it  is  important  to  know  how 
large  the  modulus  of  A  can  be  made  without  penalty. 

In  this  paper  I  present  an  extension  of  Korff  s9  approach 
for  deriving  the  transfer  function  associated  with  the  La¬ 
beyrie  technique  to  the  transfer  function  that  determines 
the  signal  of  the  KT  technique.  Neither  Korffs  result  nor 
the  general  solution  has  an  analytical  expression,  and  each 
must  be  calculated  numerically.  Results  of  such  calcula¬ 
tions  are  presented  for  a  circular,  unocculted  entrance  pupil 
of  an  aberration-free  telescope.  In  addition,  approxima¬ 
tions  for  the  cases  of  small  and  large  spatial  frequencies  are 
discussed,  which  permit  a  simpler  estimate  of  the  loss  of 
cross-spectrum  signal  as  A  is  increased. 

The  treatment  described  here  is  based  on  the  assumption 
of  log-normal  statistics  for  the  complex  wave  amplitude  per¬ 
turbations,  the  near-field  approximation,911  and  a  Kolmo¬ 
gorov  spectrum  for  atmospheric  turbulence.12  Barakat  et 
a/.13  discussed  the  influence  of  aberrations  of  atmospheric  as 
well  as  instrumental  origin  on  KT  cross-spectrum  phases  for 
a  one-dimensional  entrance  pupil.  They  assumed  a  Gauss¬ 
ian  correlation  function  to  describe  the  second-order  statis¬ 
tics  of  wave-front  perturbations.  Such  an  analysis  is  cer¬ 
tainly  adequate  for  the  study  of  the  effects  of  instrumental 
aberration.  I  largely  neglect  those  in  this  paper,  and  I  con¬ 
centrate  on  the  effects  of  atmospheric  turbulence  for  a  two- 
dimensional  entrance  pupil.  The  assumptions  for  the  wave- 
front  statistics  that  I  use  are  generally  believed  to  be  more 
relevant  to  astronomical  observations,11  and  there  is  sub¬ 
stantial  observational  evidence  to  support  them.1415  An 
excellent  treatment  of  the  statistics  of  wave-front  perturba¬ 
tions  caused  by  turbulent  media  can  be  found  in  the  review 
by  Bertolotti  et  a/.16 

2.  NOTATION 

By  using  log-normal  statistics  and  near-field  conditions,910 
the  instantaneous  optical  transfer  function  in  a  symmetrical 
formulation  may  be  expressed  as 

s-(n-/M,+? M-y) 

x  +  jj~  *,(r-  f)]}dr-  <4> 


W(r)  denotes  the  complex  amplitude  transmittance  of  the 
telescope’s  entrance  pupil  and  describes  the  telescope  aber¬ 
rations.  <t>i( r)  denotes  the  instantaneous  phase  retardation 
that  is  due  to  atmospheric  turbulence,  r  denotes  a  two- 
dimensional  pupil-plane  coordinate.  The  coordinate  f  rep¬ 
resents  a  two-dimensional  angular  frequency  with  units  of 
line  pairs  per  radian.  A  is  the  light  wavelength  and  j  =  1. 

The  statistics  of  the  random  wave-front  perturbations  are 
represented  by  the  phase  structure  function  D(p),  defined  as 

D(p)  =  ([(/>(r)  -  0(r  -  p)]2>.  (5) 

Angle  brackets  denote  ensemble  averaging.  We  shall  use 
Fried’s10  expression  for  the  structure  function,  which  is 
based  on  a  Kolmogorov  spectrum  of  atmospheric  turbu¬ 
lence12: 

D(p)  =  6.88^V/3-  (6) 

r0  is  Fried’s  seeing  parameter  and  is  related  to  the  angular 
seeing  cutoff  frequency  by  fs  =  r0/A.  Occasionally,  we  shall 
use  the  fact  that  Eq.  (6)  depends  only  on  the  modulus  |p|. 

In  order  to  simplify  the  mathematical  expressions  without 
suffering  a  loss  of  generality,  we  shall  redefine  coordinates. 
We  measure  angular  frequencies  in  units  of  the  cutoff  fre¬ 
quency  of  the  instrument,  fr  =  D/A,  where  D  is  the  entrance 
pupil  diameter,  and  we  then  set  D  =  1.  The  new  Fourier- 
space  variable  is  now  q  =  Af ID  =  Af,  and,  for  the  frequency 
region  of  interest,  0  ^  |q|  $  1,  where  Iql  represents  the 
magnitude  of  the  two-dimensional  vector  q.  Also,  Fried’s 
parameter  is  measured  in  units  of  the  telescope  diameter,  so 
we  set  a  -  r0/D.  We  can  express  Eq.  (4)  in  these  variables 
and  obtain 

S,(q)  =  j  J  W(r  +  q/2)  W*(r  -  q/2) 

X  expt/[0,(r  +  q/2)  -  0,(r  -  q/2)]|dr.  (7) 

The  domain  of  integration  is  limited  by  the  terms  W ,  which 
vanish  outside  the  area  of  the  entrance  pupil,  i.e.,  for  lr|  >  1/ 
2.  The  structure  function,  recast  in  the  new  variables,  be¬ 
comes 

D(p)  =  6.88^5/3.  (8) 


3.  DERIVATION  OF  THE  SIGNAL  TRANSFER 
FUNCTION 

The  objective  of  this  section  is  to  derive  the  ensemble  aver¬ 
age  of  the  transfer  function  CTF^(q)  of  the  cross-spectrum 
signal  by  using  Eq.  (7).  The  general  result  is  presented  in 
Subsection  3.A,  Eq.  (11).  Approximations  of  the  cases  for 
small  and  large  spatial  frequencies  are  discussed  in  Subsec¬ 
tion  3.B;  the  results  are  presented  in  Eqs.  (20)  and  (27), 
respectively.  A  description  of  the  numerical  evaluation  of 
Eq.  (11)  is  given  in  Subsection  3.C. 

A.  General  Analysis 

From  Eq.  (3),  it  is  seen  that  the  CTF  is  given  by 

CTF4(q)  =  (S,(q)S,*(q-A)>.  (9) 


'A 


m 

;$i 

v*j 

1 
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For  convenience,  we  drop  the  frame  number  index  i  in  en¬ 
semble-average  expressions.  Inserting  Eq.  (7)  into  Eq.  (9) 
yields 

CTFA(q)  =  /[  I  W(r  +  q/2)W*(r  -  q/2) 


x  i- x j *! / 1  »•  +  q/2)  -  •/>( r  —  q/2)||dr 

x  j  j  W[r'  +  (q  -  A)/2]  W(r'  -  (q  -  A)/2| 

X  expOl- 0[r'  +  (q  -  A)/2) 

+  <t>[r'  +  (q  -  A)/2]!)dr'\.  (1( 


The  general  result  for  the  transfer  function  of  KT  cross 
spectra  is  obtained  in  Appendix  A  by  extending  Korff’s9 
derivation  to  nonzero-frequency  offsets  A  The  final  result 
is  given  by 

CTFA(q)  =  |||| 

x  ^a-i-q  +  A^^ 


X  expt~V2[D(q)  +  D(q  -  A)  +  D(6  +  A/2) 

+  D(«  -  A/2)  -  D(q  -  A/2  +  5) 

-  D(q  -  A/2  -  5)])d5.  (11) 

Korffs  result  is  reproduced  in  Eq.  (11)  if  A  is  set  to  zero. 
The  first  part,  the  integral  over  a,  contains  purely  instru¬ 
mental  terms;  it  may  be  represented  as  the  overlap  area  of 
four  circles  (or  annuli)  representing  the  transparent  portions 
of  the  entrance  pupil  if  there  are  no  aberrations.  In  contrast 
to  Korffs  result,  however,  the  centers  of  the  four  pupils  are 
not  located  on  the  corners  of  a  parallelogram;  two  of  the 
pupils  are  shifted  by  the  cross-spectrum  term  A/2  (see  Fig. 
1).  Equation  (11)  has  no  analytic  solution  if  Eq.  (8)  is  used 
for  the  structure  function,  and  therefore  Eq.  (11)  must  be 
solved  numerically. 

B.  Approximations  and  Special  Cases 

1.  The  Signal  at  q  =  0 
Setting  q  =  0  in  Eq.  (11)  yields 

CTFa(0).|/|/|lv(^fS)|V(-«^) 

x  w^!±j^d<Td6  exp[—  %  D(A)]. 


If  we  restore  variables  r  and  r'  [see  Eq.  (A7)  in  Appendix  A], 
we  obtain 

CTFA(0)  =  f  f  I  W( r)l2dr  f  [  W*( r'  -  A/2) 


X  W( r'  +  A/2)dr'  exp[-  72D(A)[.  (13) 

The  first  integral  represents  merely  the  area  of  the  entrance 
pupil  if  its  transparency  is  unity.  The  second  integral  might 
be  found  to  be  a  complex  number  and  thus  might  represent  a 
systematic  error  term  if  the  telescope  has  aberrations  and 
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therefore  VF  is  complex.  If  aberrations  are  disregarded,  the 
second  integral  represents  the  overlap  area  of  two  pupils 
separated  by  A.  The  exponential  represents  the  atmospher¬ 
ic  contribution  to  the  signal  at  zero  frequency. 

To  illustrate  this  result,  let  us  assume  an  aberration-free 
telescope  with  a  circular,  unobscured  entrance  pupil,  for 
which 

'l  if  Irl  <  ‘/2 

W(r)  =  II(r)  =  •  %  if  Irl  =  \  ■  (14) 

0  if  Irl  >  % 

«. 

The  signal  at  zero  frequency  then  becomes 


CTFA(0)  =  ~  [cos-1  A  -  A(1  -  A2)1/2]exp 


Hi)  ] 


Equations  (14)  and  (15)  indicate  that,  compared  with  the 
Labeyrie  case  for  which  A  =  0,  the  signal  decreases  with 
increasing  A.  In  what  follows,  the  area-normalized  signal 
transfer  function  will  be  presented,  which  is  unity  at  the 
frequency  origin  for  the  Labeyrie  case.  This  presentation 
permits  us  to  assess  correctly  the  loss  of  signal  over  the  entire 
frequency  range  as  A  is  increased. 

2.  The  Signal  at  Small  Spatial  Frequencies 
The  regime  for  small  spatial  frequencies  Iql  «  a  can  be 
approximated  by  expanding  the  sum  of  structure  functions 
in  Eq.  (11)  into  a  power  series.  The  basic  expansion  of  the 
structure  function  [Eq.  (6)],  up  to  second  order,  is  given  by 

+  (h)!(‘-^cos!|’)  + -]}' 


Fig.  1.  Geometry  of  the  four-circle  overlap.  The  centers  of  *he 

circles  are  located  at  the  arrowheads  labeled  A,  B,  C,  and  D.  The  _ 

arrows  indicate  the  magnitude  and  direction  of  the  vectors  q,  5,  and  ,jeg 
A.  In  the  case  shown  here,  the  overlap  area  is  determined  by  only  — 
the  three  circles  A,  B,  amf-& — ^  r 
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where  6  is  the  angle  beween  the  vectors  x  and  t  and  Ixl  »  Id. 
Using  relation  (16),  the  sum  of  structure  terms  2  in  the 
exponent  in  Eq.  (11)  may  be  expressed,  up  to  first  order  in  q 
about  5  ±  A/2,  as 

V.noo  f/qVw'  /lq  — 


,  /  I S-  A/2|y/3 


Is  +  A/2 1 


mi 


8  is  the  angle  between  q  and  S  +  A/2,  and  8'  is  the  angle 
between  q  and  5  -  A/2.  In  general,  |A|  will  be  of  order  a.  If 
we  assume  that  a  «  1,  we  may  approximate  5  —  A  *  S  and  O' 
*  6  for  most  part  of  the  5  integration.  Hence,  up  to  first 
order, 

Therefore  the  turbulence  contribution  falls  out  of  the  inte¬ 
gration,  and 

CTFA(q)  »  jjjj  w(°  +  *  +  q)  W*^  +  g  ~  Q) 

x  W(^~  S~  q  +  A)  dadS 

xexpj-ad^-b^)573]}.  (19) 

The  evaluation  of  the  integrals  results  in  the  product  of  the 
aberration-limited  telescopic  transfer  function  St,a  at  spatial 
frequencies  q  and  q  -  A.  If  we  disregard  instrumental 
aberrations  and  use  Fried’s9  result  for  the  average  long- 
exposure  transfer  function  SLE(q).  we  finally  obtain 

CTFA(q)|l,l«a  *  S,,a(q)exp[-3.44^q^  Jst,0*(q  ~  A) 

*  ‘5LEfq)‘^LE^q  —  A).  (20) 

3.  The  Signal  at  Large  Spatial  Frequencies 
We  now  assume  that  ]q|  »  a,  and,  by  using  expansion  (16), 
we  find,  for  the  sum  2  of  exponentials  in  Eq.  (11)  up  to 
second  order  in  S, 

^  «  6  88  ^l«  +  A/2ljm  +  ^|S  —  A/2IJ/3  +  (qj* 


,^q~A  J'3  2^q- A/2J/3 

+  |  *-(!- i- cos’ A] 

3  |q  -  A/2|2  \  12  )\ 


+  fi - * - 15  f 1  “  *4  cos2  A  •  (21) 

3  |q  -  A/2l2  V  12  /J 

Here  8  is  the  angle  between  q  -  A/2  and  S.  Again,  A  will  be 
of  order  a  and  small  compared  with  q.  Therefore  it  is 
justified  to  assume  that  q  -  A  »  q  -  A/2  «  q.  In  this  case, 
the  terms  involving  q  cancel,  and  an  approximation  to  2,  to 
first  order  in  J,  is  given  by 


X>6.88[(l‘if/2l)“  +  (14-"21)"].  (22, 

In  order  to  proceed  further,  we  assume  that  instrumental 
aberrations  are  small  compared  with  atmospheric  wave- 
front  distortions  or,  equivalently,  that  W(r)  is  essentially 
constant  over  areas  with  a  typical  size  a.  Since  the  expon  sn- 
tial  contributes  to  the  integration  only  if  6  $  a,  the  6  varianle 
may  be  neglected  in  arguments  of  W  for  most  parts  of  the  a 
integration.  We  may  then  separate  5  integration  from  a 
integration,  and  we  obtain 

CTFa(q)  »  J  J  »"( -t  ,)lV'(”  2  ") 

+  (23, 


We  first  consider  the  a  integration  in  relation  (23).  Depend¬ 
ing  on  the  instrumental  aberrations,  this  term  may  be  com¬ 
plex  and  may  therefore  contribute  a  systematic  phase  term 
to  the  observed  cross-spectrum  phase.  If  we  neglect  instru¬ 
mental  aberrations  and  note  that  q  -  A  is  approximately 
equal  to  q,  we  obtain 

1 1 

x  w(— -q  +  A)da  <*  S0(q).  (24) 

where  So(q)  is  the  diffraction-limited  transfer  function  of 
the  instrument. 

The  integration  over  5  can  be  performed  by  using  the 
substitutions 


Fig.  2.  Decay  of  signal  in  the  cross-spectrum  transfer  as  a  function 
of  lAi/or.  Solid  line,  approximate  result  /S(A /a),  [Eq.  (26)j;  □  and  n, 
ratios  of  the  CTF  models  and  the  STF  models  lor  f  =  0  and  (  =  90 
deg,  respectively,  averaged  over  the  spatial-frequency  range  from 
—0.8  to  0.8  and  multiplied  by  0.342;  dotted-dashed  line,  measure¬ 
ment  of  C,,17  multiplied  by  0.342;  dashed  line,  measurement  of  C„17 
with  a  constant  0.1  subtracted  and  then  multiplied  by  0.342. 
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P  = 


8 

a 


n  = 


A 

2  a  ’ 


d8  =  ordp; 


(25) 


thus,  inserting  Eqs.  (25)  into  the  5  integral  and  using  Eq.  (6), 
we  obtain 


expj-3.44  ^ 


'\b  +  A/2 1 


=  «2  |  exp[— 3.44(|p  +  rj|r,/3  +  Ip  -  r,|5/3)]dp.  (26) 


The  value  of  the  integral  in  Eq.  (26)  is  a  factor  /3  that 
depends  on  A/ct.  The  numerical  calculation  of  d(A/a)  is 
presented  in  Fig.  2.  Combining  relation  (24)  and  Eq.  (26), 
we  obtain,  for  the  signal  transfer  function  in  the  high-fre¬ 
quency  approximation, 

CTF4(q)|iQ|>>rt  »  «2d^S0(q).  (27) 


4.  Combination  of  Approximate  Results 
A  simple  approximation  to  the  overall  signal  transfer  func¬ 
tion  is  given  by  combining  relations  (20)  and  (27): 

CTFA'q)  =  S,,E(q)SLE(q  -  A)  +  a2d(~jso(q).  (28) 

The  transfer  function  may  l  a  represented  by  the  sum  of  two 
components,  one  describing  the  low-frequency  portion  and 
the  other  one  describing  the  high-frequency  regime.  A  simi  - 
lar  approximation  exists  for  the  transfer  function  STF(q)  of 
the  Labeyrie  signal,11  which  is  reproduced  if  A  is  set  to  zero 
in  Eq.  (28): 

STF(q)  =  SLE2(q)  +  0.342a2S0(q).  (29) 


C.  Numerical  Calculations  of  Signal  Transfer  Functions 

Equation  (11)  was  evaluated  numerically,  assuming  an  aber¬ 
ration-free,  unocculted  circular  entrance  pupil  [Eq.  (14)]. 
The  6  integration  is  carried  out  by  a  simple  quadrature 
algorithm,  in  which  the  coordinate  6  is  represented  by  its 
magnitude  d  and  the  angle  7  that  8  makes  with  q,  so  a  polar 
coordinate  system  is  used.  The  integration  along  d  is  per¬ 
formed  first.  A  two-step  procedure  is  applied  to  determine 
the  integrand.  First,  the  value  of  the  exponential  in  Eq.  (11) 
is  calculated.  The  result  is  compared  with  the  approximate 
value  of  the  transfer  function  at  large  spatial  frequencies: 
a28(A/a).  If  the  result  is  smaller  than  a  certain  fraction  of 
the  approximate  value  (10~4  at  the  present  time)  it  is  consid¬ 
ered  insignificant,  and  the  next  step  is  omitted.  This  proce¬ 
dure  speeds  up  the  computation  considerably  while  main¬ 
taining  reasonable  accuracy  that  is  sensitive  to  the  magni¬ 
tude  of  A. 

In  the  next  step  the  a  integral  is  calculated.  The  integral 
is  represented  by  the  area  that  four  circles  of  radius  1/2  and 
with  centers  located  at  -(q  +  8)/2,  (q  -  8)/2,  (5  ~  q  +  A)/2, 
and  (8  +  q  -  A)/2  have  in  common.  The  configuration  is 
indicated  in  Fig.  1 .  Because  of  the  presence  of  the  frequency 
shift  A,  the  situation  is  more  complicated  than  in  Korffs9 
case,  and  the  symmetry  of  the  problem  is  partially  de¬ 
stroyed.  The  outer  integration  over  7  must  be  carried  out 
from  0  to  ir,  instead  of  from  0  to  ir/2  as  in  Korffs  case. 

Because  of  the  more  complicated  situation,  an  approach 
different  from  Korff s'*  was  used  in  order  to  calculate  the 
overlap  area  Impending  on  the  configuration  of  the  circle 


Fig.  3.  CTF  models  by  numerical  evaluation  of  Eq.  (11),  for  £  =  0 
deg  and  a  =  0.1.  Dotted  line,  Labeyrie  case  (STF);  solid  lines, 
models  for  (in  descending  order  of  the  curves)  IaI  =  0.02,  0.05,  0.1, 
0.15,  0.2.  For  negative  relative  wave  numbers,  the  curves  represent 
CTF  models  for  £  =  ir. 


Fig.  4.  CTF  models  by  numerical  evaluation  of  Eq.  (11),  for  |A|  = 
0.05  and  a  =  0.1.  Solid  line,  £  =  0;  dotted-dashed  line,  £  =  ?r/4; 
dashed  line,  £  =  ir/2. 


Fig.  5.  CTF  models  by  numerical  evaluation  of  Eq.  (11).  In  all 
cases,  A  =  1  and  £  =  0.  In  descending  order  of  the  curves,  «  =  0.05. 

0.1,0.15,  0.2,  0.25. 


726  .!.  Opt.  Soc.  Am.  A/Vol.  5,  No.  5/May  1988 


Oskar  von  der  Liihe 


centers,  the  overlap  area  can  be  the  area  that  two,  three,  or 
four  circles  have  in  common.  First,  the  coordinates  of  the 
circle  centers  are  calculated  for  a  given  set  of  values  for  d  and 
■).  The  two  circles  that  have  the  largest  separation  of  their 
centers  are  determined,  and  the  locations  of  their  points  of 
intersection  are  calculated.  Based  on  the  distance  of  these 
intersection  points  to  the  centers  of  the  other  two  circles,  the 
program  decides  whether  two,  three,  or  four  circles  deter¬ 
mine  the  overlap  area  and  which  ones  are  the  relevant  inter¬ 
section  points.  Finally,  the  overlap  is  computed  by  combin¬ 
ing  the  areas  of  triangles,  quadrangles,  and  circular  seg¬ 
ments.  For  the  case  of  circular  segments,  the  area  is  taken 
from  a  lookup  table  rather  than  being  computed  in  order  to 
avoid  calls  to  inverse  trigonometric  functions.  The  ap¬ 
proach  described  here  results  in  a  relatively  short  and  simple 
('  source  code  that  avoids  frequent  calculations  of  trigono¬ 
us  irn  function  values. 

Numeric  results  are  presented  in  Figs.  3-5.  Because  of 
the  frequency  shift  A,  the  transfer  functions  are  no  longer 
radially  symmetric  in  q  space,  and  they  depend  on  the  angle 
£  that  q  makes  with  A.  Figure  3  presents  the  results  for  £  = 
0;  the  curve  parameter  is  the  magnitude  of  A.  Figure  4 
presents  the  results  for  £  =  0,  for  £  =  ir/ 4,  and  for  £  =  jr/2, 
while  A  is  0.05.  In  all  cases,  a  is  0.1,  which  corresponds  to 
moderately  good  seeing  for  a  meter-class  telescope.  Figure 
5  presents  results  for  a  fixed  magnitude  of  A  and  for  varying 
seeing  parameter  a  and  demonstrates  the  loss  of  signal  with 
increasing  influence  of  turbulence. 

4.  DISCUSSION 

Figures  3-5  indicate  the  properties  of  the  signal  transfer 
function  for  KT  cross  spectra.  The  loss  of  circular  symme¬ 
try  of  the  CTF  is  caused  by  the  anisotropy  of  the  problem, 
owing  to  the  presence  of  the  frequency  shift  A;  the  CTF  has 
symmetry  of  reflection  about  the  A  direction.  The  maxi¬ 
mum  of  the  curves  moves  toward  larger  spatial  frequencies 
in  the  direction  of  A.  The  signal  decreases  dramatically  as 
lAl  increases,  under  otherwise  similar  conditions.  This  be¬ 
havior  reflects  the  decorrelation  of  the  phase  contribution  of 
the  atmosphere  to  the  instantaneous  optical  transfer  func¬ 
tion  [Eq.  (4)]  over  scales  of  order  r0/D. 

Figure  2  shows  the  loss  of  signal  as  a  function  of  | A!  in 
various  ways.  The  solid  curve  represents  the  factor  0,  which 
appears  in  the  approximation  of  the  CTF  for  large  spatial 
frequencies.  The  solid  markers  represent  suitable  averages 
of  ratios  of  model  CTF’s  for  various  values  for  A  and  the 
signal  transfer  function  calculated  for  £  =  0  and  £  =  90  deg 
and  multiplied  by  a  constant  0.342.  The  dashed  lines  show 
measurements  of  the  decrease  of  the  KT  signal  that  were 
made  by  Karo  and  Schneidermapn.17  The  data  were  taken 
from  the  curve  labeled  C,  in  Fig.  3  of  Ref.  17.  One  of  the 
dashed  curves  represents  the  data  from  Ref.  17  but  with  a 
constant  tentatively  subtracted  in  order  to  account  for  possi¬ 
ble  measurement  noise.  In  any  case,  the  numeric  models, 
the  approximation,  and  the  measurements  agree  to  a  good 
extent.  It  is  seen  from  Fig.  2  that  for  finite  A  the  signal  in 
the  KT  cross  spectra  is  always  smaller  than  the  signal  for  the 
Labeyrie  case  A  =  0.  However,  the  requirements  for  the 
magnitude  of  A  appear  to  be  less  stringent  than  was  previ¬ 
ously  believed.  In  the  approximation,  the  signal  has  half  the 
magnitude  of  the  Labeyrie  signal  if  A  =  0.6  a,  and  even  for  A 


=  a  the  signal  in  the  cross  spectra  has  some  20%  of  the 
Labeyrie  signal. 

These  results  may  be  used  to  estimate  the  influence  of 
instrumental  noise  on  the  phases  of  the  cross  spectra.  To  do 
this,  we  shall  extend  the  treatment  of  the  problem  described 
by  Deron  and  Fontanella.18  They  assumed  a  zero-mean, 
complex  Gaussian  additive-noise  contribution  with  a  vari¬ 
ance  <j2n  to  the  Fourier  transforms  of  the  individual  frames. 
They  showed  that  the  variance  of  the  cross-spectrum 
phase  4>i  of  a  single  frame  i  can  be  expressed  as 

2  <Tlro[f,(q)fl*(q-A>l 

^  ~  [lF0(q)l2STF(q)]'2  ’ 

if  the  noise-free  cross-spectrum  phase  is  arbitrarily  set  to 
zero.  Im(-)  denotes  the  imaginary  part  of  a  complex  num¬ 
ber.  It  is  assumed  in  Eq.  (30)  that  the  noise  amplitude  is 
small  compared  with  the  signal  amplitude  and  that,  because 
I  A|  is  small,  the  observed  signal  magnitude  in  the  denomina¬ 
tor  is  well  represented  by  the  signal  of  the  Labeyrie  case  A  = 
0.  I  shall  drop  the  latter  assumption  and  replace  the  denom¬ 
inator  by  the  proper  term, 


A)) 


||F0(q)F0*(q-A)|CTFA(q)]2 


(31) 


In  addition,  Deron  and  Fontanella  showed  that  the  vari¬ 
ance  of  the  imaginary  part  of  the  observed  cross  spectrum 
F,(q)F,*(q  -  A)  is  approximately 

<T?m(F,(q)F,-(q-A)|  *  «ArlF0(<l)l2STF(.j),  (32) 

and  hence 

a2  _  4lF0(q)|2STF(q) 

~  [lF0(q)F0*(q  -  A)ICTFA(q)]2 


for  a  single  frame  or 

2  ~  1  _ 4lF0(q)l2STF(q) _ 

~  N  [lF0(q)F0*(q  -  A)ICTFA(q)]2 


for  an  average  of  N  cross  spectra. 

If  we  assume  that  the  modulus  of  the  object  spectrum 
Fo(q)  does  not  vary  rapidly  over  frequency  scales  of  order  A 
(which  is  essentially  true  when  the  observed  object  has  only 
a  small  spatial  extent),  we  may  set  lF0(q)F0*(q  -  A) I  * 
|F0(q)|2,  and  we  obtain 

2  ~  _L  ^  _STF(q)_  _  (35) 

a*‘  ~  N  |F0(q)|2  [CTFA(q)]2 


A  good  representation  of  the  signal-to-noise  ratio  SNR  of 
the  Labeyrie  signal  is  given  by19 


SNR(q)  = 


|F0(q)|2STF(q) 


so  we  finally  obtain 

2  ^  _1_  1  [  STFXq) _]2 

~  N  SNR(q)  [CT FA(q)J  * 


(36) 


(37) 
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If  the  signal-to-noise  ratio  of  the  Labeyrie  signal  is  known, 
relation  (37)  may  be  used  to  estimate  the  contribution  of 
noise  to  the  error  of  cross-spectrum  phases  and  therefore  to 
the  error  in  the  reconstructed  phase.  To  illustrate  this,  let 
us  assume  a  simple,  practical  case  in  which  image  Fourier 
transforms  are  known  along  only  one  spatial-frequency  di¬ 
mension.  Data  of  this  type  are  obtained  if,  e.g.,  the  speckle 
pattern  of  a  resolved  source  in  the  focal  plane  of  an  infrared 
telescope  is  scanned  by  using  a  narrow  slit  followed  by  a 
single  infrared  detector.20  For  a  given  frequency  spacing  A, 
the  cross  spectrum  phase  </>*  is  given  at  discrete  spatial  fre¬ 
quencies  q*  =  (0,  A,  2A,  . .  . ,  kA ,  ....  K A),  where  K  is  the 
largest  integer  number  smaller  than  1/A.  Object  phases 
can  be  recovered  at  the  same  discrete  frequencies  by  apply¬ 
ing  a  simple  summation  algorithm: 

k 

tk  =  'l'k- i  +  to  = 0  or  (38) 

r  i 


(b) 

Figs.  6.  Phase  error  in  the  reconstruction  for  the  one-dimensional 
reconstruction  algorithm  [Eqs.  ( 38 ) j .  The  approximate  result  in 
Eq.  (26)  was  used  for  the  STF/CTF  ratio  in  Eq.  (37).  (a)  Phase 
error  as  a  function  of  frequency  for  IaI  =  0.02  and  a  =  0.1  and  values 
for  the  product  of  SNR  and  the  frame  number  N  of  (in  descending 
order  of  the  curves)  100, 200, 500,  and  1000.  (b)  Reconstructed  rms 
phaoe  error  averaged  overall  spatial  frtquvsci?':  .u  a  function  of  !a!, 
for  a  =  0.1.  N(SNR)  values:  A,  200; +,  500;  0,  1000 


The  variance  of  the  error  of  the  reconstructed  object  phase 
is  given  by  the  sum  of  the  error  variances  of  the  cross¬ 
spectrum  phase  values  that  were  added  to  arrive  at  a  given 
object  phase,  and  it  therefore  increases  with  increasing  fre¬ 
quency  index.  Figure  6(a)  represents  the  resulting  rms  er¬ 
ror  of  the  reconstructed  object  phase,  assuming  a  frequency- 
independent  signal-to-noise  ratio. 

The  proper  choice  of  the  magnitude  of  A  results  in  opti¬ 
mally  small  rms  phase  errors,  as  illustrated  in  Fig.  6(b).  The 
rms  object  phase  error  is  presented  as  a  function  of  A,  and  it 
is  seen  that  for  otherwise  fixed  conditions,  the  rms  error  has 
a  minimum  for  a  certain  value  of  the  frequency  shift.  For 
small  A,  the  cross-spectrum  error  is  small,  but  many  summa¬ 
tions  are  required  in  order  to  arrive  at  a  given  frequency, 
while  for  large  A,  the  cross-spectrum  errors  by  themselves 
are  large.  Although  the  optimal  choice  for  A  is  0.4a,  the  rms 
phase  error  remains  fairly  small  over  a  wide  range.  There¬ 
fore  I  conclude  that  the  magnitude  of  A  may  be  chosen 
within  the  range  from  A  =  0.2a  to  A  =  0.8a  without  introduc¬ 
ing  severe  phase  errors  due  to  random  noise  into  the  recon¬ 
struction. 

Another  application  of  the  results  presented  here  is  the 
recovery  of  the  full  object  Fourier  transform  from  the  KT 
cross  spectra.  Although  the  KT  algorithm  is  applied  mainly 
to  recover  object  phases,  object  amplitudes  can  be  obtained 
as  well  if  cross  spectra  are  calibrated  properly  for  seeing 
effects.  The  main  advantage  is  that  average  power  spectra 
no  longer  need  to  be  calculated,  which  may  be  a  matter  of 
concern  if  the  algorithm  is  to  be  applied  in  real  time.  If  a 
quick-look  reconstruction  is  all  that  is  asked  for  or  if  a 
suitable  reference  source  is  missing  but  r0  is  known,  the 
amplitude  recalibration  can  be  done  by  using  the  models 
presented  here.  The  reconstruction  algorithm  that  applies 
to  the  simple,  one-dimensional  case  discussed  above  is  given 
by 


CTFaA*-i  ’ 


K,o  =  i. 


(39) 


where  (Fi,kF’i,k-i>  represents  the  observed  cross  spectrum  at 
a  frequency  index  k,  CTFa.*  is  the  CTF  value  at  that  index, 
and  F„  k  represents  the  reconstructed  Fourier-transform  val¬ 
ue. 


APPENDIX  A:  EXTENSION  OF  KORFF’S 
ANALYSIS  TO  THE  KNOX-THOMPSON 
TECHNIQUE 

We  begin  the  analysis  with  Eq.  (10)  of  Subsection  3. A: 


CTFA(q) 


W(r  +  q/2)W*(r  -  q/2) 


X  expt/Mr  +  q/2)  -  <p( r  -  q/2)]|dr 


X 


W*[ r'  +  (q  -  A)/2]W[r'  -  (q  -  A)/2] 


X  expOl- 0[r'  +  (q  -  A)/2] 
+  <t>\ r'  +  (q  -  A)/2]|)dr'\  • 


We  rewrite  the  product  of  integrals  as  multiple  integrals, 
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regroup  factors,  and  perute  integration  and  ensemble  aver¬ 
aging: 

CTFA(q)  =  j  j  j  j  W(r  +  q/2)W*(r  —  q/2) 

X  W*[ r'  +  (q  —  A)/2l  W\r'  -  (q  -  A)/2J 

X  <expoi0(r  +  q/2)  -  4>(r  -  q/2) 

-  <t>[ r'  +  (q  -  A)/2]  +  <J>[ r‘  -  (q  -  A)/2]|))drdr'. 

(Al) 

The  next  'tep  involves  the  evaluation  of  the  ensemble 
average  in  Eq.  (Al).  Similar  expressions  were  treated  by 
Fried10  and  by  Korff.9  A  brief  justification  of  those  treat¬ 
ments  is  given  by  Roddier  (Ref.  11,  p.  293).  We  state  the 
result: 

(expO'l0(r  +  q/2)  -  0(r  -  q/2)  -  <t>[ r'  +  (q  -  A)/2] 

+  d>[r'  —  (q  —  A)/2]|)> 

=  exp(-  \  (|0(r  +  q/2)  —  <p( r  —  q/2)  —  <p[r‘  +  (q  -  A)/2] 

+  &[r'  -  (q  —  A)/2]|2)).  (A2) 

We  next  analyze  Eq.  (A2)  in  order  to  express  the  phases  <p 
in  terms  of  structure  functions.  Partial  evaluation  of  the 
square  in  the  right-hand  side  of  Eq.  (A2)  results  in 

(|d>(r  +  q/2)  -  <(>( r  -  q/2)  -  4>[r'  +  (q  -  A)/2] 

+  <t>[ r'  -  (q  -  A)/2]j2) 

=  <  Mr  +  q/2)  -  Mr  -  q/2)]2>  +  (|0[r'  +  (q  +  A)/2] 

-  0[r'  -  (q  -  A)/2] |2)  -  <20(r  +  q/2)0[r'  +  (q  -  A)/2]> 

-  (20(r  -  q/2)0[r'  +  (q  -  A)/2]> 

+  <20(r  -  q/2)0(r'  +  (q  -  A)/2)> 

+  (20(r  +  q/2)0[r' +  (q  -  A)/2]>.  (A3) 

The  first  two  terms  in  Eq.  (A3)  are  already  structure  func¬ 
tion  terms.  The  other  four  terms  may  be  rewritten  in  struc¬ 
ture  functions  if  Eq.  (A3)  is  expanded  by  the  following 
terms,  which  add  to  zero: 

( 02(r  +  q/2)  >  +  <02(r'  +  (q  -  A)/2))  +  <02(r  -  q/2)> 


CTFA(q)  = 


W(r  +  q/2)W*(r  —  q/2) 


X  W*(r'  +  (q  -  A)/2)  Wjr'  -  (q  -  A)/2] 

X  exp)—  \  [D(q)  +  D(q  -  A)  +  D(r  -  r'  +  A/2) 
+  D(r  -r  -  A/2)  -  D(q  +  r  -  r'  -  A/2) 

-  D(q  -  r  +  r' —  A/2)]|drdr'.  (A6) 


The  exponent  in  Eq.  (A6)  depends  on  only  the  difference  r  - 
r'.  A  change  of  variables  is  therefore  made: 

<j  =  r  +  r',  8  =  r  -  r', 

<r  +  6  <7-6  (A7) 

r  =  — — —  -  r  =  — ‘ 


We  then  obtain  the  final,  general  result  for  the  transfer 
function  of  KT  cross  spectra: 


CTF„<ql  -Jjjj  i|-  “) 

■K 


x  w* 


a  —  &  +  q  —  A\ 


-  6  -  q  +  A 


|  da 


2  /  V  2 

X  exp|-  \  [D(q)  +  D(q  —  A)  +  D(6  +  A/2) 
+  D(6  -  A/2)  -  D(q  -  A/2  +  6) 
-D(q-A/2-5)]|d6. 
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